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Abstract 

The three-loop form factors in massless QCD can be expressed as a linear combination 
of master integrals. Besides a number of master integrals which factorise into products of 
one-loop and two-loop integrals, one finds 16 genuine three-loop integrals. Of these, six have 
the form of a bubble insertion inside a one-loop or two-loop vertex integral. We compute all 
master integrals with these insertion topologies. 



1 Introduction 



The vertex functions of a virtual photon coupling to a quark-antiquark pair (quark form factor) 
and of a Higgs boson coupling to two gluons through an effective coupling (gluon form factor) 
are the simplest diagrams containing infrared divergencies in higher orders in massless quantum 
field theory. These form factors appear in a wide variety of applications: they can be used to 
predict the infrared pole structure of multi-leg amplitudes at a given order [1,2] and to extract 
resummation coefficients [3], and they make up the purely virtual corrections to a number of 
collider reactions (Drell-Yan process, Higgs production and decay, deep inelastic scattering). 

In the past, two- loop corrections to the massless quark [4] and gluon [5, 6] form factors 
were computed in dimensional regularisation in D = 4 — 2e dimensions to order e°. Two- loop 
corrections to this order were also obtained for massive quarks [7] . The massless two- loop form 
factors were extended to all orders in e in [8], and three-loop form factors to order e _1 (and e° for 
fermion loop contributions) were computed in [3,9]. These three- loop results had an immediate 
application in the calculation of the N 3 LO threshold-enhanced soft emission corrections [10] 
to the inclusive Drell-Yan and Higgs production cross section, demonstrating the perturbative 
stability at this order. 

In [3,9], the form factors were inferred from the behaviour of the three-loop deep inelastic 
coefficient functions [11]; this procedure can not be easily extended to yield also all the finite 
terms. Instead, one can turn to the more conventional approach of computing multi-loop Feyn- 
man amplitudes, which proceeds through a reduction [12-14] of all Feynman integrals appearing 
in the form factors to a small set of master integrals. The reduction is purely algebraic and can 
be automated using computer algebra methods [13,15]. The master integrals take the form of 
a Laurent series in e, and they must be computed to a given order in e, typically specified by 
the transcendentality of the coefficients. The finite part of the three- loop form factors requires 
transcendentality six, i.e. coefficients containing terms up to 7r 6 or £f. 

In this letter, we identify all master integrals needed for the three-loop form factors in 
Section[2j Many of these are products of integrals with one or two loops, or three-loop propagator 
integrals. Among the remaining genuine three-loop vertex integrals, several contain one-loop 
or two-loop propagator insertions. We describe how the Laurent expansion of these insertion 
topologies can be obtained either analytically or numerically in Section |3J and list the results 
for them in Section 0] Finally, Section El contains our conclusions and an outlook. 

2 Master integrals for three-loop form factors 

The topologies of the master integrals relevant to three-loop form factors can be inferred from 
two-particle cuts of the master integrals of massless four- loop off-shell propagator integrals (mass- 
less four-loop two-point functions). The master integrals of these massless four-loop two-point 
functions were identified in [16] and subsequently used in the calculation of the scalar i?-ratio [17]. 
Analytical expressions for these integrals are, however, not available in the literature. Since 
each two-particle cut is in general only one of several (two-, three-, four- and five-particle) cuts, 
knowledge of these two-point master integrals would not facilitate the calculation of the master 
integrals for the three- loop form factors. 
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These master integrals can be classified into three types: (i) products of one-loop and two- 
loop vertex functions with one off-shell and two on-shell legs, (ii) three- loop two-point functions, 
(hi) three-loop vertex functions with one off-shell and two on-shell legs. Since the one- loop and 
two- loop vertex functions are known to all orders in e [8], all master integrals of type (i) can 
be obtained directly by expansion [18, 19] of the all-orders results. Likewise, three-loop two- 
point functions appearing in type (ii) are known to sufficiently high orders in e [12,20] and are 
tabulated for example in the MINCER package [21]. The only non-trivial master integrals for 
three- loop form factors are therefore of type (hi). The full set of these integrals is displayed 
in Figure ^ Each topology contains only one master integral, which is chosen to be the scalar 
integral, with no loop momenta in the numerator and with all propagators raised to unit power. 
Nevertheless, we will give the results for the two-loop insertions for arbitrary propagator powers, 
see Section |1J The topologies A§ 2 and ^6,2 with some of the lines being massive have been 
calculated in [22] , where they enter the calculation of the three- loop matching coefficient of the 
heavy quark current. 

Among the master integrals of Figure ^ the integrals A^ t i, ^5,2, ^6,1 , ^6,3 > ^7,1 > ^7,2 are of 
special character, since they contain either a one-loop two-point insertion into a two-loop vertex 
integral (^6,3) ^7,1) ^7,2) or a two-loop (or one-loop times one-loop) two-point insertion into a 
one-loop vertex integral (^.5,1, ^5,2, ^6,1) • These so-called insertion topologies are in general 
simpler than the remaining genuine three-loop vertex integrals, since they can be obtained by 
computing a one-loop or two-loop vertex function with one or two propagators raised to a 
symbolic power. In the following, we describe the calculation of these insertion topologies. 

3 Computational methods 

Three-loop vertex integrals with one off-shell and two on-shell legs and massless propagators 
depend only on one kinematical scale: the mass q 2 of the off-shell leg. The dependence on this 
scale is given by the mass dimension of the integral, such that the coefficients of the Laurent 
expansion are constants, i.e. real numbers (which are in general of increasing transcendentality) . 
Several techniques exist to compute such single-scale integrals. 

For all one-loop and two-loop insertion topologies considered here, we performed two in- 
dependent calculations, using two different techniques: evaluation in terms of hypergeometric 
series from Feynman parametrisation and evaluation using sector decomposition. 

The Feynman parametrisation for the one-loop and two-loop vertex functions with symbolic 
powers on individual propagators results in a multiple integral in the Feynman parameters. De- 
pending on the topology, one has to integrate over at least two (one-loop vertex function) and 
at most five (non-planar two-loop vertex function) Feynman parameters. After appropriately 
decomposing the integration region to avoid parametric singularities [23], and introducing sup- 
plementary regulators at intermediate stages, one can express the results of this integration in 
terms of hypergeometric functions of unit argument, containing e-dependent coefficients. These 
can be expanded in e using the Mathematica [24] package HypExp [18] to yield the Laurent series 
of the master integrals. 

For many practical applications, and to verify the analytical results, it is sufficient to know 
the numerical values of the coefficients in the Laurent expansion of the master integrals to some 
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finite order. These can be obtained using the sector decomposition technique. 

The sector decomposition technique for the computation of multi-loop integrals is described 
in detail in [25,26]. Using this technique, the Laurent expansions of all master integrals relevant 
to the three- loop form factors can be computed to any desired order, limited only by computa- 
tion time. First applications to three-loop vertex integrals were presented already in [25]. The 
treatment of propagator powers different from unity is described in [26]. The application of 
sector decomposition to the topologies Aj t i and A-j^ has been done in two different ways: 
(a) by direct calculation of the three-loop topologies, (b) by calculating the two-loop diagram 
with e-dependent propagator powers resulting from integrating out the one-loop two-point in- 
sertion (corresponding to /5(e), Ie(e), Je(e) in Section EJ). The analytical results for A§ t i, A^^ 
and ^6,1) given for general symbolic propagator powers Vi in Section |1J also have been verified 
for some e-dependent Vi values by sector decomposition. 

The computing time for a seven propagator graph like A-? t i or A?2 U P to order e° for a 
numerical precision better than 0.1% is of the order of 20 minutes on a 2.8 GHz PC, while the 
order e term takes about 6 hours. For a precision of 1% the evaluation is about 10 times faster. 



4 Results for the insertion topologies 

In this section we list the results we obtained for the three loop master integrals with insertion 
topology. The labelling of the diagrams is according to Figure ^ The results for the diagrams 
As,i, ^5,2, and A^\ can be given for arbitrary propagator powers V{. The values of the U{ are 
assumed to be such that the arguments of all occurring T-functions are different from 

0, -1, -2,... . 

In our first diagram, namely A§ t i, we label the powers of the sloped propagators (i.e. the 
ones attached to the off-shell leg) by v\ and whereas U3, ^4, and v§ are associated with the 
three propagators that form the twofold bubble insertion. The form of the diagram immediately 
suggests that the result must be completely symmetric in \y\, v 2 } as well as in {^3, ^4, ^5}. The 
calculation leads to 



A 5 A[Vi_ 



d D k r d D i r d D 



(27r) D 7(27r) D i(27r) D P + Pi) 2 P [{k - p 2 ) 2 P [l 2 ]" 3 [(k + I + r) 2 p [r 2 ] 
l) 1 -^ r 2 . ]^/2-n r(f -i*)r(§ -Ui)T(§ -vs) 



21^5 



F (n - j£) rfa 45 - d) n*-r -n + u 1 ) r(^ - n + v 2 ) 

T(^--u U5 )r(2D-N) 



where we introduced the short-hand notations 



Vijk... = Vi + Vj + Vk + ■ ■ ■ 

N = 1^12345 • (2) 
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In the above equation (^Q), r] > is an infinitesimal quantity that indicates the way in which the 
analytical continuation has to be performed in the case q 2 > 0. 

In the special case in which all Vi are equal to unity, the result simplifies considerably. 
Defining a pre-factor Sr as 

St = Ta \ n /o t~\ / -i > ( 3 ) 



( 47r )D/ 2 



we have 



A 



5,1 Wi 



1-3 e T 6 (l - e) T(2 e) T(3 e) T(l - 3 e) 
(l-2e) (2-3e)T(3-4e) 



(4) 



In the next diagram, A^^-, the power of the upper sloped propagator is labelled by v\. v 2 
and 1/3 are the powers of the propagators of the lower bubble insertion, whereas ^4 and are 
associated with the propagators of the vertical bubble. From the form of the diagram we can 
read off that the result will be symmetric in {1^2, ^3} as well as in {^4, 1/5}. It reads 



d D k r d D i r d D 



{2v) D J {2ir) D J {2k) d P+Pi) 2 P [(l-k + P2 ) 2 r [l 2 P [{k + r) 2 r [r 2 ^ 



-1) 



l-JV 



(4^)3^/2 



3D/2-7V t(§ - v 2 ) r(f - u 3 ) r(§ - 1/4) r( 



>D 
2 



Vh ) 



r(i^)r(^)r(i^)r(i/ 4 )r(iys) 



r(jy - ffi t(d - ^45) i> 45 - g) rgf - iv + i/Q 

X F(D - u 23 )T(D - i/ 45 )r(2L> - N) 

Again, the case in which all Vi are equal to unity is much simpler, namely 

1-36 r 7 (l - e) r(e) r(3 e) T(l - 3 e) 



15,2 Wi 



1] 



-iSl 



I 7] 



(l-2e)r(2-2e)r(3-4e) 



(5) 



(6) 



The last diagram with two bubble insertions is Aq^. Again, v\ and u 2 are the powers of the 
sloped propagators. v 3 and ^4 form the powers of the upper bubble insertion, whereas 1^5 and 

are given to the lower one. The diagram also shows several symmetries, namely in {^1,^2}, 
{^3,1/4}, {^5,fg}, and, in addition, in {{^3, ^4}, {^5, vq}}- One finds 



^4-6,1 Wi] 



d D k f d u l rd u r 



(2tt) D J (2tt) D J (2tt) D l(k+Pi) 2 P [(k - p 2 ) 



21^2 



l-JV 



( 4 ^)3D/2 



i rj 



[1 2 P [(l + k) 2 p [r 2 ] U5 [(r + k) 2 ] U6 

3d/ 2 -n r(f - u 3 ) r(f - 1/4) r(f - ^ 5 ) r(f - ^) 
r( I / 1 )r(^)r(^)r( I /4)r( I / 5 )r( i / 6 ) 



t(n - gg) i> 34 - #) r> 56 - f ) r(3£ - jv+^i) r(^ _ g + V2) 

T(D-v u )T(D-v 56 )T(2D-N) 



(7) 
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where this time we have N = ^123456 • 

Finally, we again give the result for the case in which all V{ are equal to unity. 



-iSl 



1 77 



-36 T 7 (l - e) T 2 (e) T(3 e) T 2 (l - 3 e) 
r 2 (2-2e)r(2-4e) 



(8) 



Since from now on the diagrams will become more complicated, we restrain ourselves to the 
case in which the powers of all propagators are equal to unity. The remaining three diagrams to 
be considered are ^6,3, ^7 1, and each of which contains a single bubble insertion. After 
integrating out the bubble insertion we are left with an effective two-loop diagram with one 
propagator less. However, one of the propagators in the effective two- loop graph will carry a 
power that is different from unity. The two-loop crossed vertex graphs with powers different 
from unity were discussed previously in [27]. 

While computing the effective two-loop diagrams, it turns out that, after integrating over 
the loop momenta, all integrals over Feynman parameters can be carried out in a closed form. 
The respective results contain r-functions in combination with hypergeometric functions of unit 
argument. We used the aforementioned Mathematica package HypExp [18] for expanding the 
all-order results into their respective Laurent series expansions about e = 0. The explicit result 
for Aq 3 reads 



A 



6,3 



d D k f d D l f d D 



{2n) u J {2n) u J {2n) u k 2 (k-qj z (k - if (I - Pl f (r-lf r 2 

, r(Qr»(i -e) 

r(2 - 2e) 



(9) 



with 



r_iv» f d ° k f d ° l 



sl 



{2it) u J {2it) D k 2 {k-qf {k-lf (l-pi) 2 [l 2 Y 
2e T 3 (l - e) r(l - a - e) T(l - a - 2 e) 



i r\ 



r(a) r(2 - a - 2 e) T(2 - a - 3 e) 
r(l - a - 2 e) r(a + 2 e) T(a + e) T(a) T(l - a - e) 



r l-e 



r(a + 2e-l)r(l-e) 
+ ~ (1 _ 2e) -3^2(1,1-6,1-26; 2- 2e,2-a-2e; 1) 



(10) 



Substituting a = e in Eq. (|1U|) leads to the following series expansion for A 



6.3 



^6,3 = iSl 



1 rj 



-3e 



+ 



+ 



1351 55 tH 



6 

2023 
2 ~ 

26335 



6 



95 7T 



7T 935 Cs , lOvr^Cs 
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+ 



+ 



+ 65C 5 



6 6 
+1615Ca +30vr^C3 



1351 vr 2 11 7T 4 7vr 6 
+ 



18 



54 



2 2 ^3 + 585 C 5 e 3 + 0(e') 



(11) 



The above Eq. 1)10(1 can be used for two other cross-checks. First, we can consider the limit 
a — > 0. This is done by setting a = £ e, followed by the series expansion in e. Finally, we set 
£ = 0. The result has to coincide - up to a global sign - with the series expansion of the two- loop 
integral of Eq. (4) in Ref. [8]. The second check is performed by the limit a — > 1, in which 
case we have to find the result for the two-loop five propagator integral that is obtained from 
^6,3 by removing the bubble. Both checks were found to be fulfilled on the level of the series 
expansions. The calculation of /5(e) by sector decomposition provided an additional check. 
We now proceed with the integral A71, which assumes the form 



d D k r d D i r d D 



with 
I 6 (a) 



(2tt) D J (2tt) D J (2tt) D r 2 (r - k) 2 {k-qf {k-lf (k - I - p 2 ) 2 I 2 (I - P if 
r(2 - 2e) 

d D k f d D l 1 



(12) 



-(-iy /— f- 



S 2 



(2tt) u J (2tt) u [A; 2 ]" (k - q) 2 (k - l) 2 (k — I — p 2 ) I 2 (l-pi)' 
l-a-2e 



r(-e)r(2e) 



r^(l-e) r^(-e) 



4 F 3 (a,l -a-4e, l-e,-2e; 1 - 3 e, 1 - 2 e, 1 - 2 e ; 1) 



L2T(1 - 3e) 

r(l - 2 e) r(l - a - 2 e) T(2 + e) r(a + 2 e) 



+ 



r(a)r(2-e)r(l-a-4e) 
x 4 F 3 (1, 1,1 -2e,2 + e; 2,2,2 -e; 1) 
2 r(-2 e) r(l + a + 2 e) T(2 + e) T(l - a - 2 e) 

r(a)r(2-e)r(l-a-4e) 
x 4 if 3 (l, 1,1 + a + 2e,2 + e; 2,2,2 -e; 1) 

r(a + 2e)r(2-a-e) 
(l-a-2e) 2 r(a)T(2-a-3e) 
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X 4^3(1, 1 - « - 2e, 1 - a - 4e, 2 - a - e ; 2 - a - 2 e, 2 - a - 2 e, 2 - a - 3 e ; 1) 
r(-2 e) r(l + 2 e) T(2 + e) T(l + a + 2 e) T(2 - a - 2 e) 



+ ■ 



r(a) T(l - a - 4 e) T(2 - e) T(2 + 2 e) 
x 5 F 4 (l,l,2-a-2e,l + a + 2e,2 + e; 2, 2, 2 - e, 2 + 2 e ; 1) 



(13) 



Again, we have to set a = e in Eq. (|13|) in order to obtain the series expansion for A^i. It reads 

-l-3e 



A 



7.1 



i rj 



1 1 / 77 \ 1 / 77 1 

/ 27T 2 11 77 4 \ 1 

+ 4 " — — ST " 20& 7 



+ 16 



4vr 2 22 vr 4 



40C3 + 



14vr 2 C 3 



3 45 3 
8vr 2 44 vr 4 943 tt 6 



88 Cs 



45 



7560 



80 Cs 



+ 



28tt j C3 



+ 196C 3 2 - 176C 5 e + C(e 5 



(14) 



The integral Ie(a) provides another cross check since for a = 1 we have to reproduce the integral 
Aq of Eq. (5) in Ref. [8]. This we checked to be the case on the level of the series expansion. 

As we proceed, the expressions for the integrals become more and more lengthy. The result 
for the integral Aj 2 reads 



d D k f d D l fd D 



(2vr) D J (2vr) D J (2vr) D k 2 (k - qf (I - Pl f (k - I) 2 (k - I - p 2 ) 2 r 2 {r-lf 

r(e)r 3 (i- e ) 



-i Sr ■ 



■ Me) 



with 



-(-ir 



r(2 - 2e 

dPk f d D l 



(15) 



(27T) n J (27T) 1 " k [pf (k- q y (i-piY (k-iy (k-i- P2 y 

l-a-2e 



i rj 



r(l-e)r(-e)r(l-a-e) 



r(l - a - 2 e) Tja + e) T(a + 2 e) T(l - e) T(-e) T 2 (e) 
4r(a)T(l - a-4e)T(2e) 
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+ 



T(l - a-2e)r(a + e)r(a + 2e)r(-2e) 
eT(a)r(l -a-4e) 

T 2 (l - a - 2 e) r(a + e) T 2 {a + 2 e) T[a + 4e) 
r(a)T(l + 2e) 

T(l - a - 2 e) r(l - e) r(a + e) T(a + 2 e - 1) 
r(a)T(2 - a - 3e) 
x 3 F 2 (l,l - a-2e, l-a-4e;2-a-2e,2-a-3e;l) 
r(a - 1) T(l - a - 2e) T(l - e) T(-2 e) T(l + e) T(l + 2 e) 



r(a) r(l - a - 4e) T(2 - a - e) 
x 3 F 2 (1 - a,l + e, 1 + 2e; 2-a,2-a-e; 1) 
T(l - e) T(-2 e) r(l + e) T(l + 2 e) T(a + e) T(a + 2 e) 
(a + 4e) r(a) T(l - a - 3 e) T(l + a + 3 e) 
x 3 F 2 (l + e, l + 2e,a + 4e; l + a + 3e, l + a + 4e;l) 
T(l - 2e)r(l -a-2e)r(l -e)r(l + a + e)r(a + 2e) 
r(l + a)r(l-a-4e)r(2-e) 
x 4 F 3 (1, 1, 1 - 2e, 1 + a + e ; 2, 1 + a, 2 - e ; 1) 
T 2 (l -e)r(2e)r(a + e)r(a + 2e) 
(a + 2 e) r(a) T(l - a - 3 e) T(l + a + e) T(l + 2 e) 

x 4 -F 3 (l, 1, 1 - e, a + 2e ; 1-2 e, 1 + a + e, 1 + a + 2 e ; 1) 
r(l-a-2e) T 2 (l -e)r(a + 2e- 1) 
r(2-a-3e)r(l -2e) 
x 4 F 3 (l,l - a-4e,l -a-2e,l -e; 2-a-2e,2-a-3e,l-2e;l) 
a T(l - a - 2 e) T(-2 e) T(a + e) T(a + 2 e) 



+ 



+ 



+ 



+ 



+ 



r(a)T(l - a -4e) 
r(l - a - 2 e) T(-2 e) T(a + e) T(l + a + 2 e) 
r(a)r(l - a-4e) 

x 3 F 2 (l,l,l + a + 2e; 2,2; 1) 



F 2 (l,l,l + a; 2,2; 1) 



(16) 



Details about the calculation of Eq. ((Td]) can be found in Ref. [28]. Useful formulas that got 
applied at intermediate steps were taken from Refs. [29-31]. Setting a = e leads to the following 
series expansion of Aj 2 , 



A 



7,2 



i r\ 



-l-3e 



+ 4 C 3 - 



+ 



2vr 2 



+ 



83 7T 4 

360 



+ 8C3 



+ 15C5 
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/4vr 2 83^ 4 2741 vr 6 _ A 

+ 1 1 h 16 C3 

V 3 180 90720 S 



i^-73C 3 2 + 30C5]6 + O( e 2 ) 



(17) 



We finally state that the expression ()16j) for Je(a) can again be used for several cross checks. 
First, in the limit a — > 1 we have to obtain the same result as for A$ of Eq. (5) in Ref. [8] or 
of Eq. (|13|) . The check is done by first considering a = 1 + x e m d3) followed by a subsequent 
expansion in e. In the end, the limit x ~ > is carried out. A second check is provided by the 
limit a — > 0. We again set a = rje and carry out the series expansion, followed by letting r/ — > 0. 
The result has to be the same - up to a global sign - as the series expansion of is (1) of Eq. (|10|). 
All checks have been verified on the level of the respective Laurent series. 

As mentioned earlier, the coefficients of the Laurent series are real numbers. Therefore the 
method of sector decomposition is particularly well suited to compute the coefficients numeri- 
cally, thereby providing the most important check of our analytical findings. 



5 Conclusions 

In this letter, we identified and classified the master integrals required for a calculation of the 
massless three-loop quark and gluon form factors. In addition to three-loop two-point functions 
and products of one-loop and two-loop integrals, we identified 16 genuine three-loop vertex 
integrals, which are displayed in Figure ^ Among these, six integrals are so-called insertion 
graphs, containing a bubble insertion into a one-loop or two-loop vertex graph. We computed 
the master integrals for these insertion graphs analytically in a closed form which is exact to all 
orders in e, containing T-functions and hypergeometric functions. Laurent series expansions were 
subsequently obtained using the HypExp-package. All Laurent series expansions were verified 
independently using sector decomposition to determine the expansion coefficients numerically. 

The remaining ten master integrals do not contain subtopologies which would allow us to 
relate them to two-loop integrals. Their analytical computation may not be possible using 
Feynman parameters, but appears feasible with modern loop-integral techniques [32], such as 
Mellin-Barnes integration [33]. Using sector decomposition, their Laurent expansion can be 
obtained in a straightforward manner. 
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